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ABSTRACT 



It is shown that, under sufficiently intense OB-star illumination of a 
stationary photoexcitation front (PDR), nonlinear H2 photoexcitation processes 
comprising driven resonant two-photon transitions between X-state quantum 
levels, with VUV continuum light from the star supplying both driving fields, 
largely determine the photonic pathways of H2 molecules in the PDR close to 
the ionization front. Specifically, for a flux of ~4 x 10^ Habing flelds incident 
upon a PDR, the total rate at which an H2 molecule is nonlinearly photocxcited 
out of any X-state quantum level is calculated to be roughly 100 times greater 
than the total rate at which it is linearly photoexcited out of the same level. 
In strongly excited PDRs, the populations in almost all of the ~300 bound 
quantum levels of the X state will be maintained approximately equal via a 
few myriads of interconnecting two-photon steps. The remarkable importance 
of two-photon transitions in H2 photoexcitation in strongly irradiated PDRs 
derives from the exceptionally narrow Raman hnewidth (F ~10~^ sec~^) that 
characterizes all two-photon transitions between bound H2 X-state quantum 
levels. 



Subject headings: infrared: ISM ISM: clouds ISM: HII regions ISM: 

lines ISM: lines and bands ultraviolet: ISM 
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1. Introduction 

In existing theoretical models of stationary photoexcitation fronts (PDRs) [see Draine 
& Bertoldi (1996) for a comprehensive list of references], a one-dimensional geometry is 
usually considered, with light from an OB-typc illuminating star assumed incident upon 
the PDR from one direction. In all published PDR models, only linear H2 photoexcitation 
processes have been considered. The consensus among astrophysicists is that these models 
adequately describe the main photonic events that occur in H2-containing clouds separated 
by at least 10 pc from OB-type illuminating stars. In the present Letter, we draw attention 
to the fact that, in very strongly irradiated PDRs (for example, those with cloud-to-star 
distances <0.1 pc), nonlinear photoexcitation of H2 can evidently become the dominant 
photonic process determining the structure of the PDR. 



2. Resonantly Enhanced, VUV Starlight-Driven, Two-Photon Transitions 

between H2 X-State Quantum Levels 

An estimate of the resonantly-enhanced, VUV starlight-driven, two-photon transition 
rate between two H2 X-state quantum levels can be obtained with use of the formula 
shown in Fig. 1. In SI units, this formula [c.f. eq. (5.6b) of Hanna et al. (1979)] gives 
the cross-section for Stokes-wave gain in resonant stimulated Raman scattering (SRS) in a 
three-level system. However, it equally well gives the cross-section for induced nonlinear 
absorption at a;p when an intense radiation field is present at cu^ - a process often referred 
to as inverse Raman scattering (IRS). It has long been recognized (Jone & Stoicheff 1964; 
MacQuillan & Stoicheff 1966) that both processes are fundamentally the same. In this 
equation, the classical electron radius rg ~ 2.82 x 10~^^ m. The quantity Ip is the pump 
intensity in W/m^. The subscripts g,i, and / designate the ground level, intermediate 
level, and final level, respectively, of the three-level system. The quantity F is the Raman 
linewidth (in units of angular frequency), which, for the astrophysical environment assumed 
in our PDR model, would be entirely determined by a convolution of the vibrational decay 
rates of the two Xn" quantum levels involved. In light of recent calculations of quadrupole 
transition probabilities for all X-state radiative decays (Wolniewicz et al. 1998), one can 
realistically assume a value F = 10~^ sec~^ for all two-photon transitions between X-state 
quantum levels. The quantities ffi and fgi are the oscillator strengths of the Stokes-wave 
and primary frequency transitions. For all H2 Lyman (B<^=>X) and Werner (C<(=>X) 
transitions, these /- values can be simply determined from the tables in Abgrall et al. (1993 
a,b). The quantities such as Qif are the energy separations of the levels involved, expressed 
as angular frequencies, e.g. {h/2n)flif = - Ef. Note that the Stokes-wave frequency ujs 
essentially cancels Qj/ in the present case. The equation in Fig. 1 gives the cross-section in 
units of m^. 

To estimate the rates of both nonlinear and linear II2 photoexcitation in the PDR, 
one needs to know the VUV flux from the star incident upon the cloud surface. Here, 
for definiteness, we assume the illuminating star to be a BO III star, with temperature 
T=31,500K and radius R = 165^0. We take the cloud-to-star distance to be 0.1 pc, and 
assume that there is no intervening dust cloud between the star and H2-containing cloud. 
From the Planck formula for emittance of photons from a blackbody per cm^ per sec per 
unit frequency range (Allen 1976), one finds that at 1,000 A the incident flux on the PDR is 
01000 — 2.55 X 10^ photons per cm^ per sec per cm~^. Habing (1968) estimated the intensity 
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of interstellar starlight at A= 1,000 A to be Xux = 4 x 10 erg cm ^. At the surface of the 
PDR we are here considering, therefore, 0iooo is about 4.25 x 10^ Habing fields. 

As is the case in almost all published PDR models, it is here assumed that ionizing 
radiation (i.e. A < 912 A) from the star creates a thin H II region on the surface of the 
PDR cloud nearest the star. Within this H II region, a large fraction (~2/3) of the incident 
ionizing photons are converted to Lj-a photons via H"'"+ e~ recombinations followed by 
cascading transitions through H-atom excited states. Half of the Lj-a photons enter 
the neutral region, with their frequencies having become distributed in a ~20-cm"^-wide 
spectral band via the frequency redistribution that results from elastic scattering by hot 
H atoms in the H II region. One can easily estimate the total flux of ionizing photons 
incident upon the PDR from tables of the Planck function (Allen 1976). One finds this to 
be (f)ion ~ 6.2 X 10^^ photons cm~^ sec"-*^. There are thus about (f)Lya ~ 2.1 x 10^^ Ly-cc 
photons entering the neutral PDR region per cm^ per second. Thus, the Ly-a flux (per unit 

wavenumber) entering the neutral PDR region is seen to exceed that at 1,000 A by roughly 
a factor 400. In §3, we will briefly hint at an important effect this Ly-a radiation should 
have on the nonlinear excitation of H2 in strongly irradiated PDRs. 

To estimate the total rate at which an H2 molecule in the PDR would undergo 
linear photoexcitation, one can use Table 2 of Draine & Bertoldi (1996). These authors 
have calculated and summed individual contributions to the total unshielded H2 linear 
photoexcitation rate from all X^B and X^C transitions. For a flux (^1000 = 4.25 x 10^ 
Habing flelds, the total linear photoexcitation rate for a molecule in (XO, J"=l) would be 
about 1.2 X 10~^ sec~^. In the PDR model of Draine & Bertoldi (1996), a Dopplcr width 
Ril.7 cm~^ is assumed. 

With use of the equation shown in Fig. 1, we now evaluate the driven (simultaneous, 
not sequential) two-photon transition rate from XO, J"=l to XI, J"=l involving CO(n~), 
J'=l as resonant intermediate state. A rough - but conservative - approximation one can 
make in integrating the expression for ams over the applied VUV pumping continuum 
is to assign to the quantity {Qig - ujp) appearing in the denominator of the equation in 
Fig. 1 a constant value corresponding to an offset of one wavenumber, and then integrate 
over a spectral interval two wavenumbers wide on either side of the resonance. From 
the value of 0iooo earlier estimated, the value of Ip for a four-wavenumbcrs-widc spectral 
interval is found to be 0.00002 W/m^. With substitution of appropriate values for the other 
quantitities appearing in the equation of Fig. 1 (e.g. fgi = 0.03, //j = 0.07), one flnds the 
cross-section for induced absorption to be (Tirs ~ 8.4 x 10~^^ cm^. Multiplying this value 
by 01000 X 4 shows the transition rate of the IRS process to be ~8.6 x 10~^ sec~^. The rate 
of nonlinear photoexcitation of an H2 molecule in (XO, J"=l) via this particular two- photon 
transition is thus about 0.7 times the earlier estimated total rate of linear photoexcitation 
of an H2 molecule in the same (XO, J"=l) quantum level. 

However, driven two-photon transitions from (XO, J"=l) to (XI, J"=l) can also occur 
via the paired transitions Cn-OQl, Cn-lQl (n = 1-5) and also Cn-ORl, Cn-lRl (n = 0-5). 
This makes the total two-photon transition rate out of (XO, J"=l) exceed the total linear 
rate by ~7 times. (The wavelengths of both C6-0Q1 and C6-0R1 arc shorter than 912 
A, the Lyman limit. No A < 912 A starlight can penetrate beyond the HII region of the 
PDR.) One can also go from (XO, J"=l) to (XI, J"=3) via transitions Cn-ORl, Cn-lP3. 
With the inclusion of these paths, the two-photon transition rate out of (XO, J"=l) would 
be roughly 10 times the total linear rate of photoexcitation. One next must consider 
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analagous two-photon transitions from (XO, J"=l) that terminate on (X2, J"=l), (X2, 
J"=3), ...(X8, J"=l), (X8, J"=3), etc.. With inclusion of these, the two-photon excitation 
rate from (XO, J"=l) would exceed the total linear excitation rate by ~70 times. Finally, 
two-photon excitation from (XO, J"=l) can also occur via B-state resonant intermediate 
levels. Additional inclusion of such pathways, makes the total rate of two-photon excitation 
out of (XO, J"=l) exceed the total linear photoexcitation rate out of the same quantum level 
by at least a factor 100. Roughly the same nonlinear-to-linear photoexcitation rate ratio 
would obtain for H2 molecules in any of the ~300 X-state quantum levels. The estimated 
total nonlinear photoexcitation rate out of any such level is thus ~1.2 x 10~^ sec~^. 



3. Consequences of VUV Starlight-Driven Two-Photon Transitions Occurring 

in Heavily Irradiated PDRs 

From the transition rates calculated above, it is apparent that, through a few 
myriads of interconnecting two-photon transitions, the H2 populations in each of the 
~300 X-state quantum levels would be maintained approximately equal. Note that no 
H2 photodissociation can result from two-photon transitions occurring from the highest 
bound X-state levels to unbound X-state continuum levels. The two-photon transition rate 
depends inversely upon the Raman linewidth, and the latter represents a convolution of the 
lifetimes of the quantum levels connected by the two-photon process. When one steps from 
the highest bound X-state levels to unbound levels of the X-state continuum, the lifetime 
drops by roughly eighteen orders of magnitude! On the basis of this same disparity in 
lifetimes, one perhaps should question the validity of the assumption usually made in linear 
PDR models that roughly one out of ten H2 photoexcitations from X-state quantum levels 
result in photodissociation via fluorescent VUV transitions terminating on unbound levels 
of the X-state. This would be true if linear photoexcitation of H2 in PDRs did involve real 
e xcitation of B - and C-state quantum levels. However, in collisionless media, it is known 
( [Loudon 1983| ) that linear photoexcitation occurs entirely via spontaneous resonant Raman 



scattering, and, at exact resonance, the cross-section for the latter varies again as F ^ (Shen 
1974). 

If one knows the total H2 density in the neutral region of the PDR, one can 
approximately estimate the maximum thickness d of the nonlinearly excited region as 
follows. Assume an H2 density of 10^ cm~^. In the nonlinearly excited region, these 
molecules will be evenly distributed amongst all the ~300 X-state quantum levels, with the 
spontaneous infrared radiation emission rate from each level being ~10~^ sec~^. Assume 
the average IR emitted photon to have an energy of 5,000 cm~^. Equating the VUV power 
from the star incident upon the PDR in the entire spectral range of the Lyman and Werner 
bands that lies longward of 912 A to the IR power radiated in a column of thickness d: 

(20, 000)(2.55 X 10^)(100, 000) ^ (10^)(10*'^)(5, 000)(i, 

one finds d~0.003 pc. Figure 1 of [Lemaire et al. (1996)| is a striking colored image 
(colors corresponding to emission intensity) of II2 1-0 S(l) vibrational emission at 2.121 
/im in NGC 7023, taken at very high spatial resolution. It shows II2 vibrational emission 
originating from a ~0.02-pc-thick, roughly planar, cloud located ~0.07 pc to the NW of 
the illuminating star HD 200775. (From Earth, one views this cloud largely through one 
of its 'edges'.) There is very much higher H2 emission coming from three ~0.004-pc-thick 
embedded filaments located near the surface of the cloud facing the star. In view of the 



-6- 



calculations that have been presented above, we suggest that these bright filaments in NGC 
7023 might represent PDR regions that are nonhnearly excited. 

It was noted above that Lj-a radiation generated in the H II region enters the PDR 
neutral region at an intensity level roughly 400 times greater than the continuum level at 
1,000 A, but only over an estimated ~20-cm -wide bandwidth. There are a dozen or so 
allowed transitions originating from X-state quantum levels whose frequencies fall within 
such a Lj-a spectral distribution. A slight contribution towards the equalization of H2 
populations in the various X-state quantum levels arises from these resonances. However, 
there is a different - rather profound - effect that this Ly-a radiation can induce, when 
it is applied to an equally populated manifold of X-state quantum levels. It can induce 
broadband stimulated Raman scattering (SRS) to occur on the three or four strongest 
of the above mentioned resonant transitions, generating broadband IR Stokes-wave light 
on strong transitions to EF-states. An SRS process would occur as part of a 2n-wave 
parametric oscillation (SRS-PO) process, with additional IR and VUV light being generated 
on strong transitions ultimately returning molecules to the X-state levels from which the 
SRS-PO processes originated. As shown elsewhere (Sorokin & Glownia 2001), the so-called 
'unidentified infrared (emission) bands (UIBs)' can be successfully assigned to IR emissions 
coherently generated in such SRS-PO processes. 

One of the authors (PPS) acknowledges receiving strong encouragement from Dr. 
Anita J. Schell-Sorokin to publish the central idea upon which this article is based. 
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Fig. 1. — Equation for the cross-section of either SRS Stokes- wave gain or IRS- induced 
absorption in a resonant, three-level system. The light waves in the diagram are depicted 
for the SRS case. 



